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ABSTRACT 

We detect a dip of 20-45% in the surface brightness and number counts of 
NVSS sources smoothed to a few degrees at the location of the WMAP cold spot. 
The dip has structure on scales of ~ 1° to 10°. Together with independent all-sky 
wavelet analyses, our results suggest that the dip in extragalactic brightness and 
number counts and the WMAP cold spot are physically related, i.e., that the 
coincidence is neither a statistical anomaly nor a WMAP foreground correction 
problem. If the cold spot does originate from structures at modest redshifts, as 
we suggest, then there is no remaining need for non-Gaussian processes at the 
last scattering surface of the CMB to explain the cold spot. The late integrated 
Sachs-Wolfe effect, already seen statistically for NVSS source counts, can now 
be seen to operate on a single region. To create the magnitude and angular size 
of the WMAP cold spot requires a ~ 140 Mpc radius completely empty void at 
z<l along this line of sight. This is far outside the current expectations of the 
concordance cosmology, and adds to the anomalies seen in the CMB. 

Subject headings: large-scale structure of the universe - cosmic microwave back- 
ground - radio continuum: galaxies 



1. Introduction 



The detection of an extreme "cold spot" (jVielva et al. 1120041 ) in the foreground-corrected 
WMAP images was an exciting but unexpected finding. At 4° resolution, Cruz et al. (2005) 
determine an amplitude of -73 /iK, which reduces to -20 fiK at ~10° scales (ICruz et al. 
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20071 ). The non-gaussianity of this extreme region has been scrutinized, (Cayon, Jun & 
Treaster 2005; Cruz et al. 2005, 2006, 2007) concluding that it cannot be explained by either 
foreground correction problems or the normal Gaussian fluctuations of the CMB. Thus, the 
cold spot seems to require a distinct origin - either primordial or local. Across the whole sky, 
local mass tracers such as the optical Sloan Digital Sky Survey (S DSS jYork et al. 1120001 1 and 
the radio NRAO VLA Sky Surve y (NVSS. ICondon et al. I I1998T) are seen to correlate wit h 



the WMAP images of the CMB ( iPietrobon. Balbi & Marinucci 



probably through the late integrated Sachs- Wolfe effect (ISW, 



2006 



Cabre et al. 



Crittenden fc Turok 



2006) 



19961). 



McEwen et al. I ( 120071 ) extended the study of radio source/CMB correlations by per- 
forming a steerable wavelet analysis of NVSS source counts and WMAP images. They 
isolated 18 regions that, as a group, contributed a significant fraction of the total NVSS- 
ISW signal. Three of those 18 regions were additionally robust to the choice of wavelet form. 
The centroid of one of those three rob ust correlated reg ions, (#16) , is inside the 1 ° cold 
spot derived from WMAP data alone (ICruz et al. 1120071 ). although McEwen et al. I (120071 ) 
did not point out this association. 

The investiga tions reported h e re wer e conducted independently and originally without 



knowledge of the iMcEwen et al. I (120071 ) analysis. However, our work is a posteriori in 
nature, because we were specifically looking for the properties in the direction of the cold 
spot. These results thus support and quantify the NVSS properties in th e spec ific direction 



of the cold spot, but should be considered along with the lMcEwen et al. I (120071 ) analysis for 
the purposes of an unbiased proof of a WMAP association. 



2. Analysis and Characterization of the NVSS "dip" 

We examined both the number counts of NVSS sources in the direction of the WMAP 
cold spot and their smoothed brightness distribution. The NVSS 21 cm survey covers the 
sky above a declination of -40° at a resolution of 45". It has an rms noise of 0.45 mJy/beam 
and is accompanied by a catalog of sources stronger than 2.5 mJy/beam. Because of the 
short interferometric observations that went into its construction, the survey is insensitive to 
diffuse sources greater than pa 15' in extent. Convolution of the NVSS images to larger beam 
sizes, as done here, shows the integrated surface brightness of small extragalactic sources; 
this is very different than what would be observed by a single dish of equivalen t resolution. In 



the la tter case, the diffuse structure of the Milky Way Galaxy dominates (e.g, lHaslam et al 



(Il98ll )). although it is largely invisible to the interferometer. 

To explore the extragalactic radio source population in the direction of the WMAP 
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cold spot, we first show in Figure 1 the 50°x50° region around the cold spot convolved to a 
resolution of 3.4°. Here, the region of the cold spot is seen to be the faintest region on the 
image (minimum at l//,bjj = 207.8°, -56.3°). At minimum, its brightness is 14 mK below 
the mean, with an extent of «5°. The WMAP cold spot thus picks out a special region in 
the NVSS - at least within this 2500 n region. 

We examined the smoothed brightness distribution across the whole NVSS survey using 
another averaging technique that reduces the confusion from the brightest sources. We first 
pre-convolved the images to 800", which fills all the gaps between neighboring sources, and 
then calculated the median brightness in sliding boxes 3.4° on a side. The resulting image 
is shown in Figure [2] which is in galactic coordinates centered at 1/7=180°. The dark regions 
near the galactic plane are regions of the NVSS survey that were perturbed by the presence 
of very strong sources. Note that the galactic plane itself, which dominates single dish maps, 
is only detectable here between -20° < In < 90°, where there is a local increase in the number 
of small sources detectable by the interferometer. 

To evaluate the NVSS brightness properties of the cold spot, we compared it with the 
distribution of median brightnesses in two strips from this all sky map. The first strip was 
in the north, taking everything above a nominal galactic latitude of 30° (More precisely, 
we used the horizontal line in the Aitoff projection tangent to the 30° line at 1//=180°.) 
The second strip was in the south, taking everything below a nominal galactic latitude of 
-30°, but only from 10° < Iji <180°, to avoid regions near the survey limit of <5=-40°. The 
minimum brightness in the cold spot region (« 20 mK) is equal to the lowest values seen in 
the 16,800 square degree area of the two strips, and is ~ 30% below the mean (Figure [3]). 

Formally, the probability of finding this weakest NVSS spot within the pal0° (diameter) 
region of the WMAP cold s pot is 0.6%. This a pos teriori analysis thus is in agreement with 



the statistical conclusions of iMcEwen et al. I (120071 ) that the NVSS properties in this region 



are linked to those of WMAP. We also note that the magnitude of the NVSS dip is at the 
extreme, but not an outlier of the overall brightness distribution. We thus expect that less 
extreme NVSS dips would also individually correlate with WMAP cold regions, although it 
may be more difficult to separate those from the primordial fluctuations. 

The NVSS brightness dip can be seen at a number of resolutions, and there is probably 
more than one scale size present. At resolutions of 1°, 3.4°, and 10°, we find that the dip is 
~ 60%, 30% and 10% of the respective mean brightness. At 10° resolution, the NVSS deficit 
overlaps with another faint region about 10° to the west, while the average dip in brightness 
then decreases from -14 mK (at 3.4°) to -4 mK. 

The dip in NVSS brightness in the WMAP cold spot region is not due to some peculiarity 
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of the NVSS survey itself. In Figure 4, we compa re the 1° convo l ved N VSS image with the 



similar resolution, single dish 408 MHz image of lHaslam et al. I (ll98ll ). T his 408 MHz all 



sky m ap is dominated in most places by galactic emission, and was usedby iBennett et al 



( 120031 ) as a template for estimating the synchrotron contribution in CMB observations. 
On scales of 1°, the fluctuations in brightness are a combination of galactic (diffuse) and 
extragalactic (smeared small source) contributions. In the region of the cold spot, we can 
see the extragalactic contribution at 408 MHz by comparison with the smoothed NVSS 

1.4 GHz image. Note that although there is flux everywhere in the NVSS image, this is 
the "confusion" from the smoothed contribution of multiple small extragalactic sources in 
each beam, whereas the 408 MHz map has strong diffuse galactic emission as well. Strong 
brightness dips are seen in both images in the region of the WMAP cold spot - with the 
brightness dropping by as much as 62% in the smoothed NVSS; in the 408 MHz map, this 
is diluted by galactic emission. 

To look more quantitatively at the source density in the cold spot region, we measured 
the density of NVSS sources (independent of their fluxes) as a function of distance from 
the cold spot. Figure 5 shows the counts in equal area annuli around the WMAP cold spot 
down to two different flux limits. With a limit of 5 mJy, there is a 45±8 % decrease in counts 
in the 1° radius circle around the WMAP cold spot centroid. At the survey flux limit of 

2.5 mJy, the de c rease is 23±3%. This reduction in number counts is what is measured by 



McEwen et al. I ( 120071 ) in their statistical all-sky analysis. 



3. Foreground Corrections 



Several studies have claimed that the properties of the cold s pot are most likely an effect 



of in c orrect foreground su btraction (jChiang fc Naselsky 1 12004 



Coles I l2005t iLiu fc Zhang 



20051 ; iToieiro et al. 1120061). This possibility has been investigate d in detail for bot h the 



first year (IVielva et al. 1 12004] ; ICruz et al. I I2005L 120061 ) and third (ICruz et al. 1 120071 ) year 



WMAP data. The arguments against foreground subtraction errors can be summarized in 
three main points - 1) The region of the spot shows no spectral dependence in the WMAP 
data. This is consistent with the CMB and inconsistent with the known spectral behavior 
of galactic emission (as well as the SZ effect). The flat (CMB-like) spectrum is found both 
in temperature and kurtosis, as well as in real and wavelet space. 2) Foreground emission 
is found to be low in the region of the spot, making it unlikely that an over-subtraction 
could produce an apparent non-Gaussianity. 3) Similar results are found when using totall y 
independent methods to model and subtract out the foregrou nd emission (ICruz et al. 1120061) . 



namely the combined and foreground cleaned Q-V-W map (IBennett et al. 1120031 ) and the 
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weighted internal linear combination analysis (ITegmark et al. 1120031 ) 



Now that we know that there is a reduction in the extragalactic radio source contribution 
in the direction of the cold spot, we can re-examine this issue. We ask whether a 20-30% 
decrement in the local brightness of the extragalactic synchrotron emission could translate 
into a foreground subtraction problem that could generate the WMAP cold spot. We are 
not re-examining the foreground question ah initio, simply examining the plausibility that 
the deficit of NVSS sources could complicate the foreground calculations in this location. 
The characteristic brightness in the 3.4° convolved NVSS image around the cold spot is ~ 
51 niK at 1.4 GHz; the brightness of the cold spot is ~ 37 mK. This difference of 14 mK 
in brightness (4 mK at 10° resolution) represents the extragalactic population contribution 
only, as the NVSS is not sensitive to the large scale galactic synchrotron emission. By 
contrast, the single dish 1.4 GHz brightness within a few degrees of the cold spot is ~ 3.4 K, 
as measured using the Bonn Stockert 25 m telescope (IReich and Reich 1119861 1 . Therefore the 
total synchrotron contribution at 1.4 GHz is ~ 0.7 K above the CMB, 50 times larger than 
the localized extragalactic deficit. 

One way to create the cold spot would be if the universal spectral index used for the 
normal galactic (plus small extragalactic) subtraction was incorrect for the extra brightness 
temperature contribution of the NVSS dip, ST (-14 mK at L band, 1.4 GHz). We make an 
order of magnitude e stimate of this potential error. Follo wing the first year data analysis 



( IBennett et al. 1120031 ). and the similar exercise performed by lCruz et al. I (120061 ). we consider 



fitting a synchrotron template map at some reference frequency z/ re j, and extrapolating the 
model, (F(iy) mo dei), with a spectral ind ex of = — 3 to t he Q, V, and W bands. This 
spectral index is consistent with those oflCruz et al. 1 (l2006h and the averag e spectral index 



observed in the WMAP images (IBennett et al. 



2003 



Hinshaw et al. 1120061 ). Under the null 



hypothesis that the spectral index of ST is the same as that of the mean brightness T , we 
then calculate in the region of the deficit, 

F{v) mo del = [To(v n f) + S T(u re f)) (u / V re/ ) _3 . (1) 

However, if the actual spectrum of ST is -a (from L band through W band) instead of 
-3, then the true foreground subtraction, F(i/) true , should have been 



= T (u re f) {yj v re f) + ST(v ref ) {yj v ref )~ 



(2) 



The foreground subtraction would then be in error as a function of frequency as follows, 
expressed in terms of the L band temperatures : 
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6F(u) = F(u) 



true 



F(v 



model 



8T{v L ){v L /vY [1 - [v ref /vf- a ] 



(3) 



Three different r eference frequenci es have been used for synch rotron extrapolation - the 



Has lam et al. I fll981f ) 408 MHz map (Bennett et al. 1 120031). t he Ijonas. Baart fc Nicolson 



(Il998h 2326 MHz Rho des /Hart RAO survey flCruz et al. 1120061 ). and the internal K and Ka 
band WMAP images (IHinshaw et al. II2006I ). Since the foreground subtraction errors would 
be worst extrapolating from the lowest frequency template, we start at 408 MHz and look 
at the problems caused by a spectral index for ST that is different than the assumed -3. 

We obtain a rough measure of the spectral index of the dip by comparing the 1° 
resolution maps (Figure 4) at 408 MHz and 1400 MHz. At l n ,b n = 207°,-55° , we find 
5T = 2.6±0.75 K (30±12 mK) at 408 (1400) MHz, yielding a spectral index of -3.6±0.5. 
Using Equation (3), this would actually lead to a WMAP "hot spot" if a spectral index of 
-3.0 had been assumed for the extrapolation. Within the errors, the worst foreground ex- 
trapolation mistakes would then be hot spots that range from +0.5 to +4.6 fiK at Q band, 
1° resolution, (« +0.25 to +2.3 /iK at 4° resolution). Our derived spectral index for the dip 
is steeper than expected for extragalactic sources, so we also do the calculation assuming 
the flattest reasonable extragalactic spectrum of -2.5 . This would lead to a spurious cold 
spot of -2.9 /zK at 4° resolution . In either case, this is far below the -73 /iK observed at this 
resolution in WMAP, and we thus conclude that the deficit of NVSS sources does not lead 
to a significant foreground subtraction error of either sign. 



4. Discussion 



The WMAP cold spot could have three origins: a) at the last scattering surface (z 
1000), b) cosmologically local (z < 1), or c) galactic. Because the spot corresponds to a 
significant deficit of flux (and source number counts) in the NVSS, we have argued here that 
the spot is cosmologically local and hence, a localized manifestation of the late ISW effect. 



Cruz et al. I (120051 . 120071 ) derived a temperature deviation for the cold spot of ~ — 20/iK 
and a diameter of 10° using the WMAP 3 year data; on scales of 4° the average temperature is 
lower, ~ —73 /iK. Using these two data points we derive an approximate relation between the 
temperature deviation and the corresponding size of the cold spot: 6(AT/T) « 4.5 x 10~ 5 , 
where 9 is the radius of the cold spot in degrees. We now perform an order of magnitude 
calculation to see if the late ISW can produce such a spot, assuming that the entire effect 
comes from the ISW. 



The contribution of the late ISW along a given line of sight is given by (A T/T)\jsw — 



-7- 



— 2 J $ drj, where the dot represents differentiation with respect to the conformal time r/, 
drj = dt/a(t), and a is the scale factor. The integrand will be non-zero only at late times 
(z<l) when the cosmological constant becomes dynamically dominant. 

We start with the Newtonian potential given by 

$ = GM/r « -^r 2 p b 5. (4) 

The proper size r and the background density pb scale as a and a~ 3 , respectively. The growth 
of the fractional density excess, 5(a) in the linear regime is given by D(a) = 8(a)/8(a ), and 
D(a) is the linear growth factor. For redshifts below ~ 1 in ACDM, this factor can be 
approximated as 5(a) ~ a5(a )(3 — a)/2. Assuming that the region is spherical, its comoving 
radius is r c = 0.5Az(c/H), and Az is the line of sight diameter of the region. The change 
in the potential over drj can be approximated by 

A$~^r c 2 Pc 5(Az), (5) 

where subscript c refers to the average comoving size of the void and the comoving back- 
ground density. In ACDM the Hubble parameter is roughly given by H 2 = Hq(1 + 2z), 
for redshifts below ~ 1. Incorporating these approximations we get the following relation 
between the size of the region, its redshift and the temperature deviation from the late ISW: 



We now ask under what conditions this expression is consistent with the observed rela- 
tion between the size and temperature of the cold spot derived earlier, 6(AT/T) « 4.5 x 1CT 5 . 
For Q m ~ 0.3 and 8 — — 1 (i.e. a completely empty region) this leads to the simplified re- 
lation, #(1 + z) « 6, where 9 is in degrees, as before. Since the spot's association with the 
NVSS places it at z ~ 0.5 — 1, this leads to a self-consistent value of the radius of ~ 3 — 4° 
for the observed spot. For c/H = 4000 Mpc, the comoving radius of the void region is 
120-160 Mpc. 

How likely is such a large underdense region in a concordance cosmology? Suppose there 
is only one such large underdense region in the whole volume up to z=l. The correspond- 
ing void frequency is then the ratio of the comoving volume of the void to the comoving 
volume of the Universe to z=l, which is roughly 3 x 10~ 5 . Is this consistent with ACDM? 
Void statistics have been done for a number of optical galaxy surveys, as well as numerical 
structure formation simulations. Taking the most optimistic void statistics (filled dots in 
Fig. 9 of Hoyle & Vogeley, 2004) which can be approximated by logP = — (r/Mpc)/15, a 
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140 Mpc void would occur with a probability of 5 x 10~ 10 , considerably more rare than our 
estimate for our Universe (3 x 10 -5 ) based on the existence of the cold spot. One must keep 
in mind, however, that observational and numerical void probability studies are limited to 
r c ~ 30 Mpc; it is not yet clear how these should be extrapolated to r c > 100 Mpc. 

We note that llnoue & Silkl (j2006al lbh had already suggested that anomalous tempera- 
ture anisotropies in the CMB, such as the cold spot, may be explained by the ISW effect. In 
contrast to our calculation described above, their analysis considers the linear ISW plus the 
second order effects due to an expanding compensated void, partially filled with pressureless 
dust, embedded in a standard CDM (Inoue and Silk 2006a) or ACDM (Inoue and Silk 2006b) 
background. It is reassuring that the size of the void indicated by their analysis — about 200 
Mpc if located at z ~ 1 — is roughly the same as what we get here using linear ISW. 

The need for an extraordinarily large void to explain the cold spot would add to the 



l ist of a nomalies associated with the CMB. (See iHoldman. Mersini-Houghton fc Takahashi 



J2006al lbl) j'or a theory that predicts such large voids based on a particular landscape model.) 
These include the systematically higher strength of the late ISW correlat ion measured for a 



varie ty of mass tracers, compared to the WMAP predictions (see Fig. 11 of lGiannantonio et al. 



2006) , and the alignment and p lanarity of the quadrupole and the octopole (Ide Oliveira-Costa et al. 
2004uLand fc Magueijo 1120051 ) . We can, however, conclude that models linking the cold spot 
with the larger scale anomalies, such as the anisotropic Bianchi Type Vllh model of Jaffe et 
al. (2005), are no longer necessary. While we suggest that the cold spot is a loc al effect, low 
order global anisotropic models (e.g., 
needed for the low— I anomalies. 



Gumrukcuoglu. Contaldi fc Peloso 1120061 ) may still be 



5. Concluding Remarks 

We have detected a significant dip in the average surface brightness and number counts of 
radio sources from the NVSS survey at 1.4 GHz in the direction of the WMAP cold spot. The 
deficit of extragalactic sources is also seen in a single dish image at 408 MHz. Together with 
previous work, we rule out instrumental artifacts in WMAP due to foreground subtraction. 
A fuller examination of the statistical uncertainties associated with our combination of the 



McEwen et al. I (120071 ) wavelet results and our own a posteriori analysis should be performed. 
With this caveat, we conclude that the cold spot arises from effects along the line of sight, 
and not at the last scattering surface itself. Any non-gaussianity of the WMAP cold spot 
therefore would then have a local origin. 



A 140 Mpc radius, completely empty void at z<l is sufficient to create the magnitude 
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and angular size of the cold spot through the late integrated Sachs- Wolfe effect. Voids this 
large currently seem improbable in the concordance cosmology, adding to the anomalies 
associated with the CMB. 

We suggest that a closer investigation of all mass tracers would be useful to search 
for significant contributions from isolated regions. Also, if our interpretation of the cold 
spot is correct, it might be possible to detect it indirectly using Pl anck, through the lack of 



lensing-induced polarization B modes (IZaldarriaga fc Seljak 1119971 ). 
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Fig. 1.— 50° field from smoothed NVSS survey at 3.4° resolution, centered at \u, bu 
= 209°, -57°. Values range from black: 9.3 mJy/beam to white: 21.5 mJy/beam. A 10° 
diameter circle indicates the position and size of the WMAP cold spot. 
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Fig. 2. — Ait off projection of NVSS survey, centered at 1//, bu = 180°, 0°, showing the 
median brightness in shding boxes of 3.4°. The WMAP cold spot is indicated by the black 
box. Closer to the plane, large dark patches arise from sidelobes around strong NVSS sources. 
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Fig. 3. — The cumulative distribution, normalized to 1000, of median brightness levels (mK) 
in 3.4° sliding boxes of the NVSS images in two strips above |b//| > 30° (see text). The 
minimum brightness (which is from the cold spot region) is indicated by a vertical line. 
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Fig. 4. — 18° fields, with 1° resolutio n, centered at 1 



f lHaslam et al. \\198W . Right: 1.4 GHz flCondon et al 



bj 7 = 209°, -57°. Left: 408 MHz 
19981 ). A 10° diameter circle indicates 



the position and size of the WMAP cold spot. 
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Fig. 5. — Number of NVSS sources in 3.14 square degree annuli as a function of distance 
from the cold spot. The counts axis refers to the results for counts of sources with S>5mJy; 
the grey line refers to counts for S>2.5mJy with those counts multiplied by 0.56. Each bin 
is independent. 



